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In the present work we study numerically quasi-equatorial lensing by the charged, 
stationary, axially-symmetric Kerr-Sen dilaton-axion black hole in the strong deflec- 
tion limit. In this approximation we compute the magnification and the positions of 



, the relativistic images. The most outstanding effect is that the Kerr-Sen black hole 

caustics drift away from the optical axis and shift in clockwise direction with respect 
J> ■ to the Kerr caustics. The intersections of the critical curves on the equatorial plane 

o . 

as a function of the black hole angular momentum are found, and it is shown that 
they decrease with the increase of the parameter Q 2 /M. All of the lensing quantities 
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are compared to particular cases as Schwarzschild, Kerr and Gibbons-Maeda black 

O holes. 
CT 
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I. INTRODUCTION 

One of the consequences of Einstein's General Theory of Relativity is that light rays are 
deflected by gravity. Although this discovery was made in the last century, the possibility 
that there could be such a deflection had been suspected much earlier, by Newton and 
Laplace, for the first time. The phenomena resulting from the deflection of electromagnetic 
radiation in a gravitational field are referred to as gravitational lensing (GL) and an object 
causing a detectable deflection is known as a gravitational lens. The basic theory of GL was 
developed by Liebes Refsdal and Bourassa, Kantowski Detailed discussions on 
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GL are the monographs by Schneider, Ehlers, Falco 



4] and Zakharov [5j and the reviews 

and 



by Blandford and Narayan jc|, Refsdal and Surdej Q, Narayan and Bartelmann 
Wambsganss {9}. 

Nowadays, gravitational lensing is a rapidly developing area of research, and it has found 
applications ranging from the search of extrasolar planets and compact dark matter to the 
estimation of the values of the cosmological parameters id ]. In most of these cases it is 
possible to assume that the gravitational field is weak, hence the angle of deflection of the 
light in the field of a spherically symmetric body with mass M can be approximated by: 
a ~ 4GM/uc 2 , where u is the impact parameter 4J. It is well known that for a Schwarzschild 
black hole, the deflection angle diverges as u — > 3\^3M, allowing photons to wind many times 
around the lens before reaching the observer. In this case, an infinite set of images at both 
sides of the black hole appears 11]. In 1959 Darwin [12] noticed that photons passing very 
closely to a black hole may suffer much larger deflections without falling into the event 
horizon. If an observer is near the line connecting the source and the lens he would detect 
two infinite series of images very close to the black hole, besides the two classical weak field 
images. These images are produced by photons, making one or more complete loops around 
the black hole before reemerging in the observer's direction. These relativistic images could 
provide a profound test of general relativity in its strong field regime. In this situation, it 
would be surely possible to distinguish among relativistic theories of gravitation (e.g. Brans- 
Dicke, induced gravity, etc). For this reason, it is useful to perform similar investigations 
within alternative pictures and to evaluate the differences in the results. In this sense, the 
work of Bekenstein and Sanders, where they considered gravitational lensing by a cluster 
of galaxies in the frame of scalar tensor theories [13j], may be of interest. In this context 
Virbhadra, Narasimha and Chitre [l4j calculated the Einstein deflection angle for a general 
static, spherically-symmetric spacetime and constructed circular gravitational lens model 
characterized with two parameters: Schwarzschild mass and "scalar charge". The authors 
applied this results to the Einstein-Massless Scalar theory by considering the general static 
and spherically-symmetric solution, given by Janis, Newman and Winicour and proposed 
that the present lens model could be able to distinguish observationally a black hole from 
naked singularity. 

The study of the gravitational lensing when light passes very close to a massive body, 
such as a black hole (a process also known as gravitational lensing in the strong deflection 
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limit (SDL)), has received wide attention in the recent past. The development of lens- 
ing theory in the strong-field regime started with the work of Frittelli, Kling and Newman 

1 51 ] and that of Virbhadra and Ellis [16]. These authors proposed a definition of an ex- 
act lens equation without reference to background spacetime, and constructed the exact 
lens equation explicitly in the case of the Schwarzschild spacetime. They also noticed that 
the strong-field thin-lens approximation (without any small angle approximation) describes 
lensing satisfactorily even at small impact parameters. The Virbhadra-Ellis lens equation 
takes an intermediary position between the exact lens equation and the quasi-Newtonian 
approximation. It makes no assumptions to the smallness of bending angles, but it does 
make approximative assumptions to the position of the light sources and the observer. In 
order to the Virbhadra-Ellis lens equation be valid, the spacetime must be asymptotically 
flat and both the observer and the light sources must be far away from the lens. Moreover, 
Virbhadra and Ellis [16j studied lensing due to Schwarzschild black hole in an asymptotically 
flat background by numerical techniques. They defined the relativistic images which are very 
demagnified unless the observer, the lens and the source are very highly aligned and would 
serve as a proof of the general relativity if they were ever observed. After this paper Virb- 
hadra and Elli s |l7l| numerically investigated the lensing by naked singularities. On the other 
hand, Perlick 18j has considered lensing in a spherically symmetric and static spacetime, 
based on the lightlike geodesic equation without approximations. Eiroa, Romero and Torres 

19l | have described Reissner-Nordstrom black hole lensing, confirming the appearance of a 
similar patterns of images. Later on, the formulas given in Ref. were used by Petters to 
calculate relativistic effects on microlensing events. Bhadra has considered the gravitational 
lensing due to Gibbons-Maeda-Garfinkle-Horowitz-Strominger (GMGHS) charged black 
hole 211 ] . In the last two papers, the authors concentrate on the methods developed by 
Bozza. Bozza et al. {22! developed an analytical technique for obtaining the deflection angle 
in the strong-field situation and showed that the deflection angle diverges logarithmically as 
light rays approach the photon sphere of a Schwarzschild black hole. In [23!] Bozza extended 
the analytical theory of strong lensing for a general class of static spherically symmetric 
metrics and showed that the logarithmic divergence of the deflection angle at the photon 
sphere is a common feature of such spacetimes. Bozza and Mancini 24j used the strong field 
limit to study the time delay between different relativistic images, showing that different 
types of black holes are characterized by different time delay. Whisker [25] and Eiroa 26] 
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used the strong field limit approach to investigate the gravitational lensing properties of 
braneworld black holes, still Eiroa 2?J recently considered the gravitational lensing by an 
Einstein-Born-Infeld black hole. Sarkar and Bhadra have studied the strong gravitational 
lensing in the Brans-Dicke theory [28]. 

In cases of ordinary lensing, the lens is placed between the source and the observer. The 
situation when the observer is placed between the source and the lens is called retro-lensing 
and was studied for the first time by Holtz and Wheeler [29]. They analyzed only the two 
stronger images for a black hole placed in the galactic bulge with the Sun as a source and 
proposed retro-lensing as a new mechanism for searching of black holes. Eiroa and Torres 
301 ] presented a complete treatment in the strong field limit of gravitational retro-lensing 
by a static spherically symmetric compact object having a photon sphere. Their results are 
compared with those corresponding to ordinary lensing in similar strong field situations. 
Shortly after that, Zakharov at al. [3l|] discussed a possibility to detect retro-images of the 
Sun by a Schwarzschild black hole and the formation of mirages near rapidly rotating Kerr 
black hole horizons. They proposed a procedure to measure masses and rotation parameters 
through analysis of the form of the mirages. 

The SDL method was also applied to wormholes. Tejeiro and Larranaga considered 
natural wormholes and their astrophysical signatures in various cases. Applying the strong 



field limit of gravitational 



ensing theory, they calculated the deflection angle produced by 



Morris-Thorne wormholes 32j in asimptotically flat space-times and showed that wormholes 



act like convergent lenses. Recently, Nandi, Zhang and Zakharov 33] considered the SDL by 
traversable Lorentzian wormholes. Their wormhole solutions are considered in the Einstein 
minimally coupled theory and in the brane world model. The observables in both theories 
show significant differences from those arising in the Schwarzschild black hole lensing. As a 
corollary, it is shown that wormholes with zero Keplerian mass also act as ordinary lenses. 
In cosmological aspect, Mukherjee and Majumdar 34j considered particale motion and SDL 
ensing in the metric exterior of a charged dilaton black hole in de Sitter universe. Konoplya 



35| | analyzed the motion of massless and massive particles around black holes immersed 



in an asymptotically uniform magnetic field and surrounded by some mechanical structure, 
which provides the magnetic field, and showed the increases of the distance of minimal 
approach, the time delay and the bending angle for a ray of light propagating near a black 
hole. Recently Konoplya 36j considered the Ernst-Schwarzschild solution for a black hole 
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immersed in a uniform magnetic field and estimated corrections to the bending angle and 
time delay due to presence of weak magnetic fields in galaxies and between galaxies, and 
also due to influence of strong magnetic field near supermassive black holes. 

Several authors have studied gravitational lensing caused by rotating black holes. Moti- 
vation for such a study in the SDL phenomenology comes from the fact that previous works 
on the subject have selected the supermassive black hole hosted by the radio source Sgr A* 



37,|38( as the best candidate for strong field gravitational lensing as it is mentioned in Ref. 
461 ] . Unfortunately, in the first work approximations valid only for small deviations from 



the straight line path are considered, and therefore the proposed methods are not suited for 
studying relativistic images 39|, |4(|. In our aspiration to study the phenomenology of the 
relativistic images in the strong gravitational field of a Kerr black hole, no approximations 
can be taken and we need to work with the full equations of motion for null rays. There have 
been numerous articles on the motion of null rays in the gravitational field of a Kerr black 



hole 



41 



42 



43 



44] . Besides, some interesting general results have recently been derived 



through Morse theory 451 ] . 



In Ref. 



46J Bozza studied the SDL equations for quasi-equatorial orbits of photons around 



a Kerr black hole and provided analytical expressions for the positions and magnifications of 
the relativistic images. Recently after that, Vazquez and Esteban solved the lens equation 



far from the equatorial plane for some particular cases [47]. In their next work Bozza 
et. all (48] made a considerable step towards the analytical treatment of Kerr lensing, 
by restricting the observer in the equatorial plane and solving the general lens equation 
using perturbative methods for small values of the black hole spin. The limitation to the 
equatorial observer is motivated by the fact that the most important candidate for a black 
hole, Sgr A*, is likely to have a spin axis perpendicular to the galactic plane, where the solar 
system lies, in first approximation. Bozza et. all 49j generelized this work by removing 
the restriction to observers on the equatorial plane. Amore at. all in series of works 50] 



developed the analytical method for calculating the deflection angle of lig 



it in a general 



static and spherically symmetric metrics. Recently, Iyer and Petters 51] developed an 



analytical perturbation framework for calculating the bending angle of light rays lensed by a 
Schwarzschild black hole. Using a perturbation parameter given in terms of the gravitational 
radius of the black hole and the light ray's impact parameter, they determined an invariant 
series for the strong-deflection bending angle that extends beyond the standard logarithmic 
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deflection term used in the literature. 

Presently, the general relativity is not the only viable theory of gravitation. As we have 
already mentioned, the strong gravitational lensing could provide profound examination of 
ihe space-time around different kinds black holes. Therefore, following the method of Bozza 



46] in the present work we wish to study strong gravitational lensing due to a rotating, 



charge, dilaton-axion Kerr-Sen black hole 52j in the heterotic string theory with the aim 



of investigating the influence of the dilaton and the axion field in the strong deflection limit 
observations. 

The outline of this paper is as follows. The second section contains a derivation of the 
first order differential system for null geodesies applied for Kerr-Sen solution. In Sec. Ill we 
discuss the lensing equation in the equatorial plane and calculate numerically the deflection 
angle in the strong field limit. In Sec. IV we discuss the quasi-equatorial lensing and 
compute the positions of the caustic points and the magnifications of the images. In Section 
V the critical curves and the caustic structure are considered. A discussion of the results is 
given in Sec. VI. 

II. ROTATING KERR-SEN DILATON-AXION BLACK HOLE AND NULL 

GEODESICS 



In 1992 Sen [53( was able to find a charged, stationary, axially-symmetric solution |58|] of 



the field equations by using target space duality, applied to the classical Ker solution. The 



line element of this solution can be written, in generalized Boyer-Linquist coordinates 55], 
as 

, 2 / 2Mr\ , 2 2 f dr 2 I/l2 \ AMrasm 2 6 J ,, 
ds 2 = - ( 1 - —rj dt + P ( ~£ + dd ) — 2 dtd § 

( , 2 2Mra 2 sm 2 6\ . 2n7l2 
+ ir(r + r a )+a 2 + j 2 \sm 2 9d(f) 2 , (1) 

where function A and p 2 are given by 

A = r(r + r a ) - 2Mr + a 2 , (2) 

p2 _ _|_ r ^ _|_ fl 2 cog 2 q ^ 

Here M is the mass of the black hole, a is the specific angular momentum of the black hole, 
T a — Q 2 /M, Q being the electrical charge of the hole. In the particular case of a static 
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ution ([T]) coincides with the GMGHS solution investigated in 
21], while in the particular case r a = reconstructs the Kerr 



black hole, i.e. a = 0, the so 
SDL regime from Bhadra in 
solution. 

The Kerr-Sen space is characterized by a spherical event horizon, which is the biggest 
root of the equation A = and is equal to 



2M -r a + y/(2M 



4a 2 



(4) 



Hence follows that a < = 



1 



2M 2 



. Beyond this critical value of the spin there is no event 
horizon and causality violations are present in the whole spacetime, with the appearance 
of a naked singularity We will restrict to subcritical angular momenta. Radial coordinate, 
where the Killing vector is isotropic, i.e g tt = 0, define a surface called ergosphere and is 
equal to 



2M -r a + ^/(2M -r a ) 2 -Aa 2 cos 2 9 



(5) 



The null geodesic equations can be derived via the Hamilton-Jacobi equation for metric 
(1). Following work 52j, where Blaga and Blaga examined the geodesies of the Kerr-Sen 
black hole solution, we can derive the following fourth-order differential system for the null 
geodesies: 



p 2 r 



P 2 
P 2 <P 



P H 



-aE + L z suT 2 9 + -r[(r(r + r a ) + a 2 )E - aL z 



A 

E(r(r + r«) + a 2 ) 2 - 2MraL z 
A 



- a 2 E sin 2 9, 



(6) 
(7) 
(8) 

(9) 



where 



R(r) = [aL z -(r(r + r a )+a 2 )E) 2 -A[(L z -aE) 2 + IC], 
9(0) = K- cot 2 9[L 2 Z - sin 2 9E 2 a 2 }. 



(10) 
(11) 



In Eqs. OH])- (Q, the dot indicates the derivative with respect to some affine parameter, in 
Eqs. (JTU1) and ( fTTJ) /C is a separation constant of motion, E = —pt is the energy at infinity, 
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L z = is the angular momentum with respect to the rotation axis of the black hole, while 
a = L z /M is the angular momentum per unit mass. 

The final expression of the lightlike geodesies follows from Eqs. (EI)- ([9]) in terms of inte- 
grals 

dr f e d6 

' (12) 



±^W) J ±VWY 

A0 = a I Mr + r a)+ pE-aL a ] f [L M ^0-_oE\ 

1 ±a^/W) J ±VW) 



The sign of ^jR(r) and a/©(#) is positive when the lower integration limit is smaller 
than the upper limit, and negative otherwise. 

We describe the light ray trajectory with the assumption that at infinity, where there is 
no gravitational field, it is a straight line. In this treatment we can identify the approximate 
light ray with three parameters. The first of them, the inclination ip , is the angle that the 
incoming light ray forms with the equatorial plane, defined by 9 = ir/2. The projection of 
the trajectory in the equatorial plane has an impact parameter u. 

The distance between that point of the projection which is closer to the black hole and 
the trajectory itself we will denote by h. For the observer with coordinates (r D , o ) in the 
Boyer-Lindquist system we can define two celestial coordinates £i, £2 for the images. The 
coordinate £1 describes the observable distance of the image with respect to the symmetry 
axis in direction normal to ray of sight, while the coordinate £2 represents the observable 
distance from the image to the source projection in the equatorial plane in the direction 
orthogonal to the ray of sight. In this formulation using Eqs. ([6]), ([7j) and (|HJ) and further 
allowing E = 1 and r D — ► 00 we get 

£1 = r^sin6> : 



dr 
2 d6 



Jsm^Oo, (14) 
hsm6 . (15) 



Taking under consideration the fact that 9 Q = n/2 — ip Q and setting £j = u, it is possible 
to relate the constants of motion J and K, to the initial conditions characterizing the light 
ray: 

J = ucosipc (16) 
JC = h 2 cos 2 ip + (u 2 — a 2 ) sin 2 ip Q . (17) 
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III. LENSING IN THE EQUATORIAL PLANE 

A. Equation of the gravitational lens in the equatorial plane 

Let us consider the case when both the observer and the source lie in the equatorial 
plane of the Kerr-Sen black hole and the whole trajectory of the photon is limited on the 
same plane. If we set the black hole in the origin, then the angle between the direction of 
the source and the optical axis will be denoted by 7. 7 ~ is the case of almost perfect 



alignment. As it is well known (see Ref. 46]) the caustic points can be expressed with the 
strong field limit coefficients and can be estimated numerically as a function of the black 
hole spin. From the lensing geometry we can write the relation: 

7 = -a(Q) + 6 + 6 mod27r, (18) 

where according to this, when u <C (D LS , Dql) one can write 

J^LS ^LS 

where 9 is the impact angle of the source, Dls is the lens-source distance and Dol is the 
observer-lens distance. 

The lens equation in the equatorial plane is then 

7 = D ol + Dls _ a ^ mod ^ (2Q) 
Dls 

The angle 7 can take any value in the interval [— tt, 7r]. a(9) is the deflection angle on the 
looping light ray and it depends directly (through the impact parameter u) on the curvature 
of the space-time geometry of the rotating black hole. As we shall see below, after calculation 
of the a(8) for the Kerr-Sen black hole in the strong field limit we can solve the equatorial 
lens equation in order to derive the relativistic images. 

B. Deflection angle in the equatorial plane 

Let us start with the conditions 9 = tt/2 and h = ip = 0, which set the light ray on the 
equatorial plane. In this case, if we substitute x = r/2M to a new radial coordinate and 
measure all distances in units 2M = 1 we obtain the reduced metric in the form 

ds 2 = -A(x)dt 2 + B(x)dx 2 + C(x)d(j) 2 - D(x)dtd(f). (21) 
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The metric coefficients are: 



A(x) 
B(x) 

C(x) 
D(x) 



30 | X q 



X \~ X 

x(x + x a ) 
x(x + x a ) — x + a 2 

x(x + x a ) + a 2 + 



X \~ X ry 



2a 



(22) 
(23) 

(24) 
(25) 



X I X (y 

where x a = Q 2 /2M 2 . 

We express the derivative t and </> in terms of the metric coefficients A(x), C(x) and D(x 
Using the equations for the constants of motions E = —p t = —Qt^ and J 
written for metric (21), we obtain 

" D + 2AJ ' 



9^ 



AAC + D 2 
2C — DJ 
AAC + D 2 



(26) 
(27) 



Moreover, the Lagrangian £=^g fJiU x^x u is another constant of motion, which vanishes for 
null geodesies and can be used to express x: 



x 



±2i 



'C- J(D + AJ) 



(28) 



B(AAC + D 2 ) ' 

Since in the case under consideration h — tp — 0, immediately from Eq. ( fTBl) follows that 
the angular momentum J is equal to the impact parameter u. Evaluating the Lagrangian 
at the minimum distance of approach of the light ray, where x — 0, with the help of Eqs. 
(|26|) . (|27|) we find the precise relation between J and the closest approach distance x 

-D + y/4A C + D 2 



J = u 



2Ar 



-a + (x + x a )^/xl + x (x a - 1) + a 2 



(29) 
(30) 



Xq + X a — 1 

where all metric function are evaluated at x = xq. The sign before the square root is chosen 
to be positive when the light ray is winding counterclockwise. For a > the black hole 
rotates counterclockwise, while for a < the black hole and the photons rotate in converse 
direction. 
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FIG. 1: The radius of the photon sphere as a function of the Kerr-Sen black hole angular momenta 
for x a equal to (black, solid line), 0.02, 0.08, 0.18, 0.32, 0.5, 0.72, 0.98 (wine, short, dashed, 
dotted line). 



A starting point of the strong field limit expansion is the photon sphere, which has been 
denned by Virbhadra and Ellis Ji| and has subsequently been investigated by Claudel, 
Virbhadra and Ellis 56(. Combining the null geodesies equations and Eqs. ( I26i) -( l29l) we 
obtain the photon sphere equation for reduced stationary, axially-symmetric metric (12ip 



A C - A' C + J(A' D - A D' ) = 0, 



(31) 



for which we require to admit at least one positive solution. The real root external to the 
horizon of this equation defines the radius of the photon sphere x ps = r ps /2M. For the 
Kerr-Sen metric the photon sphere equation takes the form 



Ax\ + I2{x a - l)x 3 + {13x 2 a ~ 22x Q + 9)x 2 + 2(3x 3 a - 6x 2 a + 3x a - Aa 2 )x 



+ %l - 2X1 + X a ~ 4 « ^ 



(32) 



The radius of the photon sphere is computed numerically and is plotted in Fig. 1 for different 
values of x a . For all photons, a critical value of the angular momenta a cr exists at which the 
photons fall on the ergosphere x es = r es /2M. In the equatorial case, for the radius of the 
ergosphere we have When the photon sphere coincides with the ergosphere 

(i.e x ps = 1 — x a ), we can calculate the critical value a cr via the equation of the photon 
sphere. The values of a cr for different 
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2 Q 

&cr 



0.02 0.08 0.18 0.32 0.50 0.72 0.98 
0.354 0.348 0.332 0.304 0.262 0.204 0.124 0.01 



When a — > | 1 — x a \ then x ps — > xh, i-e the photon sphere coincides with the horizon in 
the limit of extremal Kerr-Sen black hole. 

Having the expression (26) and (28) we find the azimuthal shift as a function of the 
distance 

d± = ± ^B\Aq\(D + 2JA) (33) 

dx yfCy/AAC + D^sgn(A )[A -A% + ±(AD - A D)\ ' 

where the fact that matric coefficient A is negative under the ergosphere is taken under 
consideration. 

From the symmetry in the phase of approach and departure we can write the whole 
deflection angle via integration of Eq. (31) from Xq to infinity 

a(x ) = (j>f(x ) - it, (34) 
/ (xo) = 2^°°^x. (35) 

(f)f(xo) is the total azimuthal angle. With the decreases of the distance of closest approach 
Xq the deflection angle increases, and for a certain value of Xq the deflection angle becomes 
2tt, so that the light makes a complete loop around the black hole. Let Xq decreases further, 
then the light ray will wind several times around the lens before reaching the observer and 
finally when xq becomes equal to the radius of the photon sphere x ps the deflection angle will 
become unboundedly large and the photon will be captured by the lens object. Besides, if 
the photons are winding in the same direction of rotation as the black hole, a much smaller 
photon spheres exists for them in comparison to the photon spheres of conversely winding 
photons. The retrograde photons may be captured more easily from the black hole than the 
direct photons. 

We can find the behaviour of the deflection angle very close to the photon sphere following 



the evaluation technique for the integral (1351) developed by Bozza 23] . The divergent integral 
is first splitted into two parts one of which <f)^(xo) contains the divergence and the other 
4>^(xq) is regular. Both pieces are expanded around Xq = x ps and with sufficiently large 
accuracy are approximated with the leading terms. At first, we express the integrand of 
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( |35|) as a function of two new variables y and z which are defined by 

y = A(x), (36) 

z = . 37 

l-l/o 

where yo = A . 

The whole azimuthal angle than takes the form 

<f ) f( x o) = R{z,x )f(z,x )dz, (38) 
J o 

where the functions are defined as follows 

f(z,x ) = 1 (40) 

^sgn{A )[Ao - A% + j(AD - A ,D)] 

All functions in (139]) and fj40|) without the subscript "0" are evaluated at x = A~ 1 [(l — y )z + 
yo]. The function R(z,xq) is regular for all values of z and xq but /(z,xo) diverges when 
z — > 0, i. e. as one approaches the photon sphere. The integral (1351) is then separated in 
two parts 

<f> f {x ) = <pf{x ) + <f)f{x ), (41) 

where 

<f>f{x ) = [ R(0,x ps )fo(z,x o ) (42) 



contains the divergence and 



<f>f(xo) = / g(z,x )dz (43) 



is a regular integral in z and Xo- To find the order of divergence of the integrand, we expand 
the argument of the square root of f(z, x ) to second order in z and get the function fo(z, %o)'- 

fo(z,x ) = — 1 , (44) 
\J az + pz 1 



where: 



a = sgn{A Q y l ^\ a q C' q - A' C + J(A' D - A D' )}, (45) 
^o^o 

P = 8 MAo) { \ J 2 ^ 2 {MoCoA* + (C C' ' - 2C' 2 )y A' - C C' y A" 

+J[A C (A'^D' - AqDq) + 2A' C (A Q D' - A' D )]}. (46) 
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The function g(z,Xo) is the difference between the original integrand and the divergent 
integrand 

g(z,x ) = R(z,x Q )f(z,x ) - R(0,x p8 )fo(z,x ). (47) 

When xq becomes equal to x ps the equation of photon sphere holds. Then, the coefficient a 
vanishes and the leading term of the divergence in f is z~ x . Therefore the integral diverges 
logarithmically. The coefficient (3 takes the form 

(1 — A s) 2 , , „ n 

(3 ps = s 9 n (A ps )—— — 2 [A ps C ps — Ap S C ps + J(A ps D ps — A ps D ps )}. (48) 

ZLy ps /± ps 



Close to the divergence, Bozza [23[ obtained the analytical expression of the deflection an- 
gle expanding the two parts of the original integral (38) around Xo = x ps and approximating 
the leading terms. The result is 

a(9) = -am (?E°L _ l) + l + 0(u-u ps ), (49) 
V u ps j 

where the coefficients a, b and u ps depend on the metric function evaluated numerically at 
x ps . Dol is the distance between the lens and the observer, and 6 = u/Dol is the angular 
separation between the lens and the image. The strong deflection limit coefficients of the 
expansion (f49|) are 

-D ps + yj4A ps Cp S + D^ s 
u ps — TTk > l DU J 

^■ft-ps 

R(0,x ps ) 

a— — , (51) 

2 v/ ftps 

1 u ps I A ps | ( D ps + 2 J A ps ) J* 
where bn is the regular integral (fif (x ) evaluated at the point follows 



1 



g(z,x ps )dz. (53) 



For the same metrics the coefficient b^ cannot be obtained analytically. In such cases it 
must be evaluated numerically. 

Fig. 2 shows the strong deflection limit coefficients as functions of a. The minimum 
impact parameter has similar behaviour as the photon sphere x ps and decreases with a. a 



15 




FIG. 2: Strong deflection limit coefficients as a function of the Kerr-Sen black hole angular mo- 
menta for x a equal to (black, solid line), 0.02, 0.08, 0.18, 0.32, 0.5, 0.72, 0.98 (wine, short, dashed, 
dotted line). 

grows, while b decreases for all values on x a . Both coefficients diverge when a — > | 1 — x a \. 
The divergence of the coefficients of the expansion means that the bending angle in the 
strong deflection limit fHUl) no longer represents a reliable description in the regime of high 
a. 
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C. Relativistic images 

Once the deflection angle is known, the positions of the images can be obtained from Eq. 
( I20p . We consider a generic geometric disposition when the source, the lens and the observer 
are not aligned. 

Let us start with solving the lens equation when 9 = u/Dql ^ 1- Then, from the reduced 
lens equation 7 = —a(9®) mod 2ir and the expression for the deflection angle derived in the 
strong field limit (149]) . the positions of the n-th relativistic images can be found in first order 



46] 



0° n = ^(l + e n ), (54) 

b+~i-2Trn /~~\ 

e n = e a , (55) 

where n=l, 2, ... indicates the number of winding done by the light ray around the black 
hole. To build a correction of this solution we can expand ot(9) around 9® 



i _^OL^ e _ e oj mod2vr . (56) 



6°n 

aDr 



Replacing this result in the lensing equation ( 1201) and neglecting the hig her order terms, we 



find with sufficient precision the position of n-th relativistic images 



46] 



9 n ^8:(l- U ^° L T t DLs) ). (57) 

According to the past oriented light ray which starts from the observer and finishes at 
the source the resulting images stand on the eastern side of the black hole for direct photons 
(a > 0) and are described by positive 7. Retrograde photons (a < 0) have images on the 
western side of the black hole and are described by negative values of 7. 



IV. QUASI-EQUATORIAL LENSING 



In order to obtain the caustic structure and to calculate the magnification of the images, 
following work 46], we need a two dimensional lens equation. Following this plan, we shall 
restrict for simplicity our consideration to small values of the initial data ipo, and we shall 
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first analyze the ray trajectory close to the equatorial plane for the Kerr-Sen black hole, 
following the approach of Bozza for quasi-equatorial motion in the strong deflection limit 
mentioned above. For small declination we have a precession of the photon orbit. 



A. Precession of the orbit for small declinations 



We describe the light trajectory with an additional coordinate that is equivalent to the 
declination ip = | — 9. For small declinations and small heights h compared to the projected 
impact parameter u we have ?/>o ~ h/u. From Eqs. (TI6"1) - (117I) with the accuracy of the first 
order terms we get 

J ~ u, (58) 
JC~h 2 + u 2 ip 2 , (59) 
u = \Ju 2 — a 2 . (60) 

During the motion of the photons we get a simple evolution equation for ip as a function 



of the azimuthal angle for the Kerr-Sen black hole 46] 

^ = (6i) 

where 



— h 2 

V> = y=2+^o 2 , (62) 

_ _ a 2 + x(4>)(x(4>) +x a - 1) 
U[(t>) - U [a + J(x^)+x a -l)}x^y [b6) 

The solution of Eq. flBTl) . since x depends on <f>, is 

V(0) =^cos(0 + o ), (64) 

with 

r <t> 

u{(j)')d(j)' . (65) 



For the photons coming from infinity and returning to infinity according to strong deflec- 
tion limit the upper integral can be written as follows 



/ R w (z,x )f(z,x Q )dz, (66) 
Jo 
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FIG. 3: The coefficients a and 6 as a function of a Kerr-Sen black hole angular momentum for 
x a equal to (black, solid line), 0.02, 0.08, 0.18, 0.32, 0.5, 0.72, 0.98 (wine, short, dashed, dotted 
line). 



where 



R u (z,x ) = lu(x)R(z,x ) 



(67) 



is a regular function. The functions R(z,Xq) and f(z,Xo) are given respectively via Eqs. 
(EH) and (HO). 

Then the representation of <fif in the strong field limit is 



u 



ps 



ps 



b = b p + a In 



4:f3pg Cps 



iipslApsl {Dps ~\~ 2t/^4p S ) 



where 



i 

bR / [Ruj^Z y Xp s *) f (Z , Xp s ) - R LU (0,x ps )fo(z,x ps )}dz. 
'0 



(68) 
(69) 
(70) 

(71) 



The quantities a and 6 are evaluated numerically in the photon sphere x ps . For the 
Kerr-Sen metric the coefficient a can be calculated analytically and the result is 



1 + 



x ps [ X ps (p^ps 1) 2a(a £) -\- x a \3Xp S -\- 3(x a l)xp S -|- (x a 1 



(72) 
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where 

f = \Jx ps {x ps + x a -l) + a 2 . (73) 

In the particular case x a = 0, the Kerr solution is recovered because of which the coeffi- 
cient a is equal to 1 for all values of the black hole spin. The behaviour of a and b for all 
possible values of the black hole spin is plotted in Fig. 3. The coefficient a has a finite value 
in the extremal black hole limit , while b diverges. 

The analysis of the integrand function u(<p) leads to u < 1 for a > and to uo > 1 for 
a < 0. In the particular case of Schwarzschild black hole when a — x a — 0, u — 1. In the 
case when u < 1, (pj < <ftf and after each loop the direct photons lag with A(p towards the 
initial declination if}. When u > 1, 0j > 0/ and the retro photons reach the same initial 
declination as advances with A<f> before completing a loop. For Schwarzschild black hole a 
precession is not present. 



B. Lensing at small declinations 



When we want to study the caustic structures and to compute the magnification of the 
images for small displacement from the equatorial plane we can write the polar lens equation 
in addition to the equatorial lens equation. In order to consider the lens geometry within 
the small declination hypothesis, we can set the source at height h$ on the equatorial plane. 
We are setting the observer on height ho and during the whole time we will assume that 
the following relations 

u < (h , h 8 ) < (Dol, D LS ) (74) 
hold. After neglecting the higher order terms, using (l?4"j) we can build the polar lens equation, 



which connects the positions of the source and the observer 46] 



h. 



h 



o 



■ sm . 



u 



7 



COS ( 



7 



fn n N 1 D lD ls . - 
(Dol + D LS ) cos(p f sm . 



u 



(75) 



Our purpose is to determine the inclination if> under which the observer sees the light 
ray. The solution of (1751) for the n-th image is 

(h s + h cos (f> f n )u- h Q D LS sin0 / n 



D lD ls sin (f) f>n - u(D OL + D LS ) cos < 



(76) 



7,n- 
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The quantity (f)j n is the phase of the n-th image and may be obtained from Eq. fl68|) once 
we know the solution 9 n of the equatorial lens equation. 

According to the lensing geometry, in the polar direction we can write the parameter h 
of the incoming light trajectory as follows 

h = h Q + D OL %- (77) 

In the approximation f Cl and h <m the height h n of the n-th image is 

h = u{h s D OL - h D LS cos0 /;n ) ^ 

DolD L s sin / n - u(D l + D LS ) cos (f) fjU 

Following ( 1741) we can see that in the neighborhood of <pj = kn the denominators of the 
■0o,n an d h n vanish and both quantities diverge. Then equation 

K(j) = u{D OL + D LS ) cos 4> f - DolD L s sin 4> f = (79) 

defines the position of the caustic points. From Eq. (|68l) for the phase 0j and the formula 
for the deflection angle (T491 and with the help of the equatorial lens equation 7 = — a(#°) 
mod 2ir the caustic points equation for Kerr-Sen black hole takes the form 

+ • ' kn. (80) 



The solution 

Ik = -b+*(b-kn) (81) 
a 

shows that for each k one caustic point for direct photons and one caustic point for retrograde 
photons exists. For k = 1 we have weak field caustic points and an azimuthal shift about ir 
and strong field limit caustic point when k > 2. 

What is the magnification of the created images formed by the source close to the caustic 
points? The magnification, defined as the ratio of the angular area element of the image 
and the corresponding angular element of the source that the observer would see if there 
were no lens, in the considered lens geometry is 

d 2 Ax {D OL + D LS f 1 

" = d?As = DTs \J\- (82) 

From the lens map, which has the form 

7 = 7(0), (83) 
h S = MMo), (84) 
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FIG. 4: The angular position of the first five relativistic caustic points: k = 2, 3, 4, 5, 6 from below 
to above as a function of the Kerr-Sen black hole angular momentum for x a equal to (black, 
solid lines), 0.08 , 0.32, 0.72 (purple, short, dotted lines). 



we can write the Jacobian determinant 
c?7 dhs 



J 



dedipo 



with 



aD OL \u(D OL + D LS ) cos<p f - D olD ls sin <j> f \ 



6+7 



(85) 



(86) 



where 7 codes the number of loops made by the photon and take negative values, besides 
7 mod 27T is the angular position of the source with respect to the lens position, n = 
is the number of loops done by the photon. In the particular case 7 = — 2m7r, the source 
is situated behind the lens, while, when 7 = —(2m + l)ir, the source is situated before the 
lens. 

In Fig. H] the positions of the first five relativistic caustic points are plotted as a function 
of the Kerr-Sen black hole angular momentum for different values of x Q . All caustic points 
coincide in the situation a = x a = when the Schwarzschild metric is recovered. In this 
case 7^ = ~{k — Then, the first relativistic, caustic point of the source before the lens is 
72 = —it and as follow 73 = —2tt is the caustic point of the source behind the lens. All caustic 
points are anticipated for negative a and delayed for positive a. When x a grows the caustic 
curves become ever more shifted from one to other, while at large angular momentum they 
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move very far from the initial position. In the vicinity of an extremal black hole angular 
momentum, where the strong deflection limit approximation fails the relativistic caustic 
points diverge. Therefore, the graphics are truncated in the vicinity of the separable point 
where the divergence starts. 

From Eq. (1851) we can see that the Jacobi determinant diverges in the caustic points. 
Therefore, in order to write the magnification for the enhanced images created by a source 
close to the caustic points we can expand Eq. f!79|) around 7^ 



K{i) ~ K'( 7fc )(7 - 7.) = - hD ° L _ DLS { 1 - lk {a)). (87) 

(X 

The formula for the magnification of the enhanced images then takes the form 

"* " D 2 OL D 2 LS | 7 -7,|' (88) 

u(a)u ps (a)e lk{a) 
^k\ a ) -2 > l°yj 

CI 

b — klT 

e 7fe (a) = e a . (90) 

The quantity Jl k in essence represents the magnifying power of the Kerr-Sen black hole close 
to caustic points. 

In Fig. 0the numerically evaluated magnifying powers /7 2 to Jl % at the caustics points from 
72 to 76 respectively are plotted. The magnification decreases with growth of \a\ linearly for 
smaller values of x a , only grows for bigger plotted values of x a and diverges for all plotted 
values of x a when a approaches the extremal value 1 1 1 — x a \. The magnifications of the 
neighbouring enhanced images falls rapidly with the winding of the photons increasing ever 
more with x a and diverges at the extremal value of the black holes angular momentum, 
where the strong deflection limit approximation fails. 

The shapes of Jl k are identical for every k, but with factor (190]) they are distinguished as 
follows 

~ e -f. (91) 
For a Kerr-Sen black hole the ratio depends on x a and the angular momentum, and is 



i-'k 



plotted in Fig. [6] 
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FIG. 5: The magnifying power at the caustic point 72, 73, 74, 75, 76 as a function of Kerr-Sen 
black hole angular momentum for x a equal to (black, solid line), 0.08 , 0.32, 0.5, 0.72 (purple, 
short, dotted line) from below to above. 
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FIG. 6: The ratio of the magnifying power at a neighbouring caustic points as a function of 
the Kerr-Sen black hole angular momentum for x a equal to (black, solid line), 0.08 , 0.32, 0.72 
(purple, short, dot line) from below to above. 

V. CRITICAL CURVES AND CAUSTIC STRUCTURE 



In the basic case of Schwarzschild black hole when the light source, the lens and the 
observer are strictly aligned we have a large weak field Einstein ring and an infinite series 
of concentric relativistic Einstein rings, very close to the minimum impact angle 9^ 23]. 

For the Kerr-Sen metric in the quasi-equatorial approximation in the strong deflection 
limit we can calculate the points of the intersections of the critical curve with the equatorial 
plane, which are 

(D L + D LS ) 



ot r - o°k cr 



ciDolDls 



9 



0,cr 



ps 



D, 



l + e. 



OL 



Ik, 



(92) 
(93) 



nole in the center of our Galaxy 
38j point M BH = 3.6 x 1O 6 M 



The numerical results are plotted in the Fig. [7] for the black 
according to the current estimates for the black hole mass 
and assuming that Dol = 8.5 kpc and Dls = 1 kpc. For every value of x a we have a 
family of characteristics plotted for k from 2 to 6, where for k from 3 to 6 the curves become 
very close to each other. The critical curves diverge at the extremal value of the angular 
momentum, when the strong deflection limit approximation breaks down. Close to the point 
of the divergence the graphics are not shown. For all values of x a the critical points are 
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FIG. 7: Intersections of the strong deflection limit critical curves with the equatorial plane for 
k from 2 to 6 (from above to below) as a function of a Kerr-Sen black hole angular momentum 
for x a equal to (black, solid lines), 0.18, 0.32, 0.5, 0.72 (purple, short, dotted lines). With an 
accuracy of the order of 0.01 fi arc sec the critical curve points #jr r coincide with the positions of 
the innermost images 6^ for k > 4. 

closer to the optical axis for positive a in case of direct photons and further from the optical 
axis for the negative a when we consider opposite photons. For the static case when the 
GMGHS black hole is realized the critical points are not shifted. Therefore, as seen from 
the side of the observer, the critical curves are distorted and shifted towards the western 
side, if north is the direction of the spin. With an accuracy of the order of 0.01 fi arc sec 
the critical curve points 6^ coincide with the position of the innermost images #oo for k > 4. 
From this point of view Fig. [7] also shows the positions of the outermost relativistic images 

#00 ■ 

7fc(— \a\) and Tfe(l a l) represents the intersections of the fc-th caustics with the equatorial 
plane. For the static case (a = 0) of GMGHS black hole the caustic points coincide with the 
caustic points for the Schwardzschild black hole and they are distributed over the optical 
axis as follows 7^ = —{k — l)ir. For a ^ the relativistic caustics 72, 74, 76 are shifted 
toward the western side, while the relativistic caustic points 73, 75 are shifted towards the 
eastern side. The projections of the caustic curves on the equatorial plane are plotted in 
Fig. [HI [H [10], as seen from the north direction in order to disentangle the influence of the 
two parameters x a and a over the caustics. Generally the behaviour of the caustics as a 
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FIG. 8: The first six caustics of the Kerr-Sen black hole lens for a = 0.1, marked by arcs between 
7/t(— \a\) and 7fc(|a|) for x a equal to (black, solid lines), 0.02, 0.08, 0.18, 0.32, 0.5 (navy, short, 
dashed lines). 




O 



FIG. 9: The first six caustics of the Kerr-Sen black hole lens for x a = 0, marked by arcs between 
7fc(— \a\) and 7fc(|a|) for angular momentum a equal to 0.05 (black, solid lines), 0.09, 0.13, 0.17 
(black, short, dashed, dotted lines). 



function of the black hole angular momentum and the parameter x a are the same. For all 
values of x a and the black hole angular momentum a the relativistic caustics drift in the 
clockwise direction, while the non-relativistic caustics 71 stay close to the optical axis only 
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FIG. 10: The first six caustics of the Kerr-Sen black hole lens for x a = 0.08, marked by arcs 
between 7fc(— \a\) and 7fc(|a|) for angular momentum a equal to 0.05 (orange, solid lines), 0.09, 
0.13, 0.17 (orange, short, dashed, dotted lines). 

for small value of x a and a, and drift in counterclockwise direction with the increase of x a 
and a. With the increase of k the caustics become larger and further in comparison to their 
initial position with respect to optical axis. Moreover, as x a and a grow, the caustics drift 
in clockwise direction with respect to the caustics for smaller x a and a. From the numerical 
results in the Fig. H]we can see that at height angular momentum and x a the caustics may 
become enormous. There the strong deflection limit fails and we have not considered this 
situation. 



VI. DISCUSSIONS AND SUMMARY 

In this paper we discuss the features of light propagation close to the equatorial plane of 
the rotating dilaton-axion black hole. Modelling the massive compact object at the center 
of the galaxy as a Kerr-Sen black hole, we estimate the numerical values of different strong- 
lensing parameters. When compared to the corresponding lensing observables due to the 
Kerr black hole, it is found that there is a significant dilaton-axion effect present on the 
observable parameters. 

As the Kerr-Sen solution of the field equation describes a charged, stationary, axially- 
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symmetric black hole we can obtain all lensing phenomena for the GMGHS charged, 
static spherically-symmetric black hole when the angular momentum a = and for the 
Schwarzschild black hole when a = x a = 0. A generic source not aligned to the optical axis 
gives rise to extremely faint relativistic images, while a point source aligned to the opti- 
cal axis produces infinitely bright images. For static-spherically symmetric black holes the 
caustics are points which are disposed over the optical axis, therefore a source close to one 
caustic is simultaneously close to all other caustics and gives rise only to enhanced images. 

In the case a ^ and x a = 0, when the Kerr lensing is recovered, the caustics drift away 
from the optical axis and only one image of all other composed infinite series of images can be 
enhanced, while the rest of the images will be quite faint. If we set the source on the caustic 
point 7fc(|o|) then the outermost relativistic image on the western side will be enhanced, 
while if the source is at 7&(— then the enhanced image will be the first on the eastern 
side. For the source disposition at the caustic points 7^+1 (|a|) or jk+i(~ M) the enhanced 
relativistic image will be the second image respectively on the western or the eastern side. 
If the source is inside the caustic, the parity of the corresponding image changes parity 
and two additional images appear with the same parity as the original image. The created 
images move rapidly in vertical direction and they are not liable to description through the 
quasi-equatorial approximation [23l. 157]. 

In Kerr-Sen lensing, when we have a 7^ 0, x a 7^ the situation reminds of the Kerr 
lensing with several correction. First of all, the equatorial critical points 6^ decrease with 
the increase of x a for the whole spectrum of a. The second consequence is that a correction 
of the position of relativistic caustic points exist in the case x a 7^ 0, which influences the 
caustic lines 7&, as they become longer with the increases of x a and black hole angular 
momentum a. The third consequence is that the magnifying power for the fc-th caustic 
point, becomes larger with the increase of x a for k > 2, unlike the case k = 2, where the 
magnifying power decreases with the increase of x a . Moreover for k = 2 (the first strong 
caustic), the magnifying power decreases with the increase of a for every x a , while from 
k = 3 to k = 6, a value x' a exists at which j2k does not depend on the angular momentum. 
For all values x a > x' a , the magnifying power only grows. The fourth consequence is that, 
the ratio of the magnifying power at neighbouring caustic points depends on the angular 
momentum when x a 7^ and grows with x a as the function . 

The caustics become longer with the increases of x a and/or the black hole momentum a. 
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They also drift away from the optical axis and shift in clockwise direction in comparison to 
the Kerr caustics. Through a = the source is simultaneously close to all caustics and gives 
rise only to the enhanced images, while for the Kerr-Sen black hole the source can be close 
to only one caustic at the same time and thus produces only one enhanced image. Photons 
which move in the equatorial plane in direction opposite to the direction of rotation of the 
black hole have a bigger impact parameter compared to the impact parameter of the direct 
photons, which means that the retrograde photons are captured more easily by the black 
hole than the direct photons. The probability for absorbtion for the two kinds of photons 
decreases when the parameter x a grows. The magnification of the image strongly depends 
on the charge of the black hole and increases with x a . 

When comparing these results to the corresponding observable lensing quantities due to 
the Schwarzschild and Kerr black holes, it is found that there is a significant dilaton-axion 
field effect present in the observable parameters. This means that the mentioned black holes 
solutions are quite different with respect to the observational effects in terms of the strong 
deflection limit. Hence, detecting the relativistic images, which might be possible in the 
near future through the very long baseline interferometry high resolution imaging, we will 
be able to either accommodate or rule out the Kerr-Sen black hole candidate as a possible 
lensing object. 
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